Carrier dynamics in epitaxial graphene close to the Dirac point 
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We study the carrier dynamics in epitaxially grown graphene in the range of photon energies 
from 10 - 250 meV. The experiments complemented by microscopic modeling reveal that the carrier 
relaxation is significantly slowed down as the photon energy is tuned to values below the optical 
phonon frequency, however, owing to the presence of hot carriers, optical phonon emission is still the 
predominant relaxation process. For photon energies about twice the value of the Fermi energy, a 
transition from pump-induced transmission to pump-induced absorption occurs due to the interplay 
of interband and intraband processes. 

PACS numbers: 78.67. Wj, 63.22.Rc, 78.47.jd 
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Graphene, consisting of a single atomic layer of car- 
bon atoms in a hexagonal lattice, exhibits a unique band 
structure with zero energy gap and linear energy disper- 
sion, the so-called Dirac cone. The band structure gives 
rise to several remarkable properties, some of which are 
highly attractive for novel optoelectronic and photonic 
devices [l[ . Of key importance are the dynamics of elec- 
tronic relaxation, in particular due to the interaction with 
the phonon system. During the last years many insights 
into the relaxation dynamics have been obtained from 
single-color and two-color [6l-fTojj pump-probe exper- 
iments. Thermalization of the nonequilibrium electron 
distribution via electron-electron scattering on a sub- 100 
fs timcscale and efficient scattering via optical phonons 
on a 100 fs - few ps timescale have been identified. Com- 
mon to previous pump-probe experiments is an excita- 
tion energy of ~1.5 cV, i.e. high above the Dirac point. 
Also most graphene photonic applications demonstrated 
so far involve near infrared or visible light [ll, 12 1. The 
nature of graphene as a gapless material with constant 
absorption in the range, where the band structure is well 
described by Dirac cones, suggests to expand the studies 
into the mid infrared and terahertz range. In particular, 
it is crucial to obtain thorough insights into the relax- 
ation dynamics in the range of the optical phonon energy 
and below, i.e., close to the Dirac point, where Coulomb 
as well as optical and acoustic phonon processes can be 
significant and were both interband and intraband pro- 
cesses are relevant. Here graphene serves as a model 
system to understand the relevance of electron-electron 
and elcctron-phonon interaction for both intraband and 
interband relaxation in materials with small or vanishing 
energy gap. 

In this Letter, we study the carrier dynamics close to 
the Dirac point by varying the excitation energy E by 
more than an order of magnitude (245 meV - 10 meV). 
An optical phonon bottleneck is observed in the range E 



= 245 - 30 meV, with decay times increasing from sub-ps 
to several 100 ps. Microscopic calculations based on the 
density matrix formalism explain these time constants by 
optical phonon scattering and additionally reveal contri- 
butions due to Coulomb- and acoustic phonon-induccd 
processes. For E < 30 meV, a striking and unexpected 
sign reversal of the pump-probe signal is found. The in- 
volved processes are understood by a simple model based 
on the dynamic conductivity including intra- and inter- 
band absorption. 

We study multilayer graphene containing ^70 layers 
grown by thermal decom pos ition on the carbon termi- 
nated surface of 4H-SiC [13j. It has been shown that 
such samples behave graphcnc-likc rather than graphitic 
due to rotational stacking of layers which preserves the 
electronic symmetry of a single graphene layer 14|. The 
graphene-like nature of our samples has been verified 
by Raman spectroscopy [l5[ and magneto-spectroscopy 
3, l3 • In contrast to a recent two-color experiment on 
similar samples, where parts of the signal could be selec- 
tively attributed to the few highly doped graphene lay- 
ers at the interface our single-color study with much 
lower excitation energies investigates the relaxation dy- 
namics in the large number of almost intrinsic graphene 
layers on top of the interface layers. Since our photon 
energies are below twice the value of the Fermi energy 
of the interface layers, no interband transitions in these 
layers are excited. The free-electron laser FELBE serves 
as a source for tunable picosecond pulses of infrared radi- 
ation [73 ■ FELBE operates in a wavelength range from 
4 — 280 fim (photon energy E = 310 - 4 meV) with a rep- 
etition rate of 13 MHz. In the experiments (performed 
at E = 245 meV, 72 meV, 51 meV, 30 meV, 20 meV, 
18 meV, 14 meV, 8 meV [19() the pump beam was me- 
chanically chopped, the polarization of the probe beam 
was rotated by 90° with respect to the pump beam and a 
polarization analyzer was placed in front of the detector 
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FIG. 1: (Color online) Experiment: Normalized pump- 
induced transmission for different photon energies E (a) mea- 
sured at 10 K. Illustration of the optical pumping (arrows 
in colors like in part (a)) in comparison with relaxation via 
optical phonons (b). Temperature dependence of the pump- 
induced transmission for E = 245 meV (c) and E = 72 meV 
(d) . The symbols are experimental data, the lines are fits tak- 
ing into account the laser pulse duration and a bi- exponential 
decay. 



TABLE I: Full width at half maximum (FWHM) of the FEL 
pulse and the decay components of the pump-probe signals 
measured at T = 10 K. For small photon energies, n is below 
the experimental time resolution. 



E(meV) 


Tpulsc (ps) Tl (ps) 


T 2 (ps) 


73 (ps) 


245 


0.7 0.5 ±0.1 


5.2 ±0.2 


within noise 


72 


3 


14.5 ± 1 


300 ± 50 


51 


11 


25 ±2 


300 ± 50 


30 


7 


25 ± 2 


300 ± 50 



in order to minimize scattered pump radiation. Mercury- 
cadmium-telluride detectors optimized for the respective 
wavelength were used for the three highest photon ener- 
gies. For the lower photon energies a liquid He cooled 
extrinsic Ge:Ga detector was employed. 

In Fig. 1(a) the pump- induced transmission change 
obtained for different photon energies down to 30 meV 
is shown. The lattice temperature was 10 K. The dif- 
ferent rise times correspond to different durations of the 
FELBE pulses (see Tab. 1). The secondary peaks ob- 
served in the signals for the two highest photon energies 
are artifacts due to a second pump pulse which is caused 
by internal reflections in beamline windows used for the 
respective wavelengths. To extract relaxation times the 
experimental data were fitted taking into account a Gaus- 
sian excitation pulse and a bi-exponential decay. The re- 
sulting time constants are summarized in Tab. 1: For E 



= 245 meV, we observe a fast (0.5 ps) and a slower (5.2 
ps) relaxation time. The observation of two timescales is 
consistent with previous results obtained in experimen- 
tal studies in the near infrared showing a fast component 
ranging from 70 - 120 fs and a slower phonon- induced 
component in the range of 0.4 - 2.5 ps [2j-|4[. Our room 
temperature value T2 = 3.4 ± 1 ps obtained for E = 245 
meV is only slightly longer, thus pointing to an optical- 
phonon related mechanism. Note that at this excitation 
energy the relaxation is due to efficient interband scat- 
tering on optical phonons of an energy of -~200 meV be- 
cause direct intraband relaxation is strongly suppressed 
(cf. Fig. 1(b)). In contrast, for near- infrared excita- 
tion, intraband relaxation is the dominating relaxation 
channel. As the excitation energy is reduced to values 
below the optical phonon energy, ti is strongly increased 
to 14 - 25 ps. Additionally an even longer timescale of 
T3 = 300 ± 50 ps is observed, which was buried in the 
noise floor for E = 245 meV. Furthermore, we find an 
interesting temperature dependence: While for E = 245 
meV both the maximum transmission change and the re- 
laxation times depend only weakly on temperature, both 
quantities drop sharply with increasing temperature for 
E = 72 meV (cf. Fig. 1 (c)-(d)). At room temperature 
the relaxation for E = 72 meV is equally fast as for E = 
245 meV within the measurement accuracy. 

To understand the involved relaxation d yna mics more 
thoroughly, we use a microscopic theory [20J to calcu- 
late the relaxation dynamics of non-equilibrium carriers 
resolved in time, momentum, and angle. Within the den- 
sity matrix formalism, we derive the graphene Bloch- 
equations [20M22I ] describing the momentum dependent 
carrier density pAk (in the valence and conduction band) , 
the microscopic polarization, and the phonon occupation 
number taking into account both optical and accoustic 
phonons at the T and K points. The Coulomb and the 
carricr-phonon interactions are considered using second 
order Born-Markov approximation resulting in a micro- 
scopic Boltzmann-like scattering equation with time and 
momentum dependent scattering rates and explicitly in- 
cluding Pauli blocking terms. To compare our micro- 
scopic results with the experiment, we approximate the 
differential transmission signal by the change of the car- 
rier density in the optically excited state Apk (23|. In 
accord with the experiment we find that the relaxation is 
significantly slowed down as the photon energy is reduced 
below the optical phonon energy (Fig. 2(a)). Also the 
calculated temperature dependence (Fig.2(b)-(c)) con- 
firms the experimental observations (Fig.l(c)-(d)). How- 
ever, in contrast straight-forward expectations [24| . op- 
tical phonons are still the dominant process for cooling 
the electronic system, even if the excitation energy is 
below the optical phonon energy. This is illustrated in 
Fig. 2(c), where the curve for T = 300 K is also calcu- 
lated with the optical-phonon contribution switched off. 
The reason for this is the strong carrier-carrier interac- 
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FIG. 2: (Color online) Theory: Normalized change of the 
carrier density of optically excited states for different photon 
energies at f0 K (a). Temperature dependence of the relax- 
ation dynamics at E =245 meV (b) and E — 72 meV (c) for 
T = 100 K, 200 K, 300 K. The dashed line is calculated with- 
out optical phonon modes illustrating their importance even 
at the low energy. 
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FIG. 3: (Color online) Experimental pump-induced transmis- 
sion change for three different fluences for the photon energy 
30 meV (a) and 20 meV (b). Calculated pump-induced trans- 
mission of a single graphene layer for photon energies of 30 
meV (c) and 20 meV (d) obtained with the simple model 
based on electron heating described in the text. 



tion in graphene which provides a hot carrier distribu- 
tion on sub-ps timcscales, allowing hot electrons to scat- 
ter efficiently with optical phonons. This is energetically 
possible because some fraction of the hot electrons have 
energies exceeding the the optical-phonon energy. The 
scattering by acoustic phonons is slower and leads to a 
decay time T3 w 300 ps, which is also consistent with the 
experiment. For E = 245 meV the calculation reveals 
another feature seen in the experiment, namely an initial 
fast time constant t± . The microscopic calculation allows 
us to identify this contribution as Auger-type processes: 
the predominant impact excitation redistributes the en- 
ergy of the excitation pulse to an increasing number of 
charge carriers leading to a fast cooling of the electronic 
system. Unlike the initial thermalization process, where 
only conduction band electrons are involved, the impact 



ionization transfers electrons form the valence band into 
the conduction band [2(| . While the microscopic calcula- 
tions give an excellent qualitative agreement with exper- 
imental results, the quantitative values at low tempera- 
ture for T\ and T2 are by about a factor of 2 larger in 
the theory, which could be due to neglected non-Markov 
channels or scattering via defects. The experimentally 
observed effect of increased temperature on both the re- 
laxation time and the maximum transmission change is 
stronger (especially for E = 72 meV) than in the cal- 
culation. One possible reason for this is the neglected 
momentum dependence in the electron-phonon coupling 
matrix elements and the phonon dispersion, which is here 
considered linear for acoustic phonons and constant for 
optical phonons. Furthermore the role of the SiC sub- 
strate is only included via a global dielectric screening 
constant. The effects discussed here are expected to re- 
sult in modified relaxation times especially at elevated 
temperatures, where they become more prominent due 
to the broader thermal distribution. 

Now we turn to discussing the carrier dynamics in the 
vicinity of 2Ef of the almost intrinsic graphene layers. 
As E is decreased from 30 to 20 meV, a striking change 
in the sign of the pump-probe signals is observed (Fig. 3 
(a) and (b). The positive pump-probe signal, i.e. pump- 
induced transmission, is caused by the bleaching of the 
interband transition (cf. right inset of Fig. 4) similar to 
the cases discussed above. If the photon energy becomes 
smaller than twice the value of the Fermi energy, this 
process is not possible anymore. However, the carrier 
distribution can be heated via free-carrier absorption. As 
a result, states for which interband absorption is possible 
get populated with electrons (or equivalenty, depleted of 
holes) as indicated in the left inset of Fig. 4 25j. Hence, 
pumping in this regime results in induced absorption. For 
a quantitative description, we consider the real part of 
the dynamic conductivity [ToL l2fj| 



BT fln(2cosh(~E f /2k B T el )) 



2 



irh(ujT 
tanh 



/hu + 2E f 
V 4fc B T e; 



tanh 



fhu- 2E f 
V 4fc B T e; 



(1) 



Here <jq = e 2 /Ah is the "universal value" of the dynamic 
conductivity, T e i the carrier temperature and r the mo- 
mentum relaxation time. The first summand of eq. (1) 
indicates intraband Drude absorption, the second one in- 
terband absorption. In Fig. 4 the dynamic conductivity 
is plotted for Ef = -13 meV and t = 300 fs. This value for 
the momentum relaxation time has been deduced from 



the linewidth in magneto-spectroscopy experiments [17 1 . 
Even at 10 K the Drude peak at low energies merges with 
the tail of the interband absorption, i.e. there is finite ab- 
sorption at all energies. The arrows indicate the increase 
(decrease) of the absorption with increasing temperature 
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FIG. 4: (Color online) The dynamic conductivity for different 
temperatures (black: 10 K, blue: 20 K, green: 50 K and red 
100 K). The inset sketches the carrier distribution before and 
after the pump beam for E < 2\Ef\ and E > 2\Ef\. 



for the photon energy 20 meV (30 meV). For comparison 
with the experimental data of Fig. 3(a)-(b), we con- 
sider heating and relaxation of the electron distribution 
via the equation dT e i/dt = P(t)/c p - (T e i - T; Qtt )/T 2 , 
where P(t) = I(t)Re(a(T e i))/eoc is the absorbed power 
per area (J, laser intensity; c, speed of light), c p the spe- 
cific heat [28[ , and r 2 = 25 ps a phenomenological energy 
relaxation time (see Tab.l). This simple model cannot 
account for the full dynamics; however, it contains the 
most important processes, since, at the timcscales rele- 
vant here, a thcrmalizcd carrier distribution can be ex- 
pected. The calculated transmission changes for a sin- 
gle graphene layer are plotted in Fig. 3 (c) and (d). 
The model captures the experimentally observed sign re- 
versal. The agreement with the experiment is best for 
\Ef | = 13 ± 2 meV. From magneto-spectroscopy experi- 
ments of samples grown with the same method \Ef\ = 8 
meV has been concluded [l^l ■ We note that the interplay 
of interband and intraband processes has been suggested 
recently to explain a change in sign of the pump-probe 
signal as pump fluence is varied in a UV-pump near- 
infrarcd-probc experiment 27 1. 

In summary, we have studied the relaxation dynamics 
of carriers in graphene close to the Dirac point. A slow- 
ing down of relaxation has been observed for energies 
below the optical phonon energy, but optical phonons 
still play an important role, as shown through micro- 
scopic calculations. The theoretical description of the 
results shows that the interplay of electron-electron and 
electron-phonon scattering produce new and efficient re- 
laxation channels in situations when both intraband and 
interband relaxation are allowed, as they are in many 
narrow or vanishing gap systems. Most remarkably, the 
sign of the pump-probe signal is inverted, as the photon 
energy is varied across the value of 2Ef. This behavior 
can be understood by the interplay between inter- and 
intraband absorption and, besides being a fascinating be- 
havior by itself, it allows a rather accurate determination 



of the Fermi energy. Furthermore the effect can be the 
basis of innovative THz optical switches and active filters. 
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